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Abstract 

We show that 1-loop quantum corrections to the potential energy 
density in supcrsymmctric hybrid inflation, outside the inflationary 
valley, can not be neglected. A method is presented to calculate these 
1-loop corrections and they are applied to the case of D-term hybrid 
inflation, where a significant amount of inflation is shown to occur after 
spontaneous symmetry breaking. Taking this into account improves 
the agreement with WMAP measurements. A gauge coupling of up 
to 0.3 is still consistent with the CMB density perturbation. The 
spectral index is predicted in between 0.98 and 1.00 and the cosmic 
string contribution to the CMB anisotropy is sufficiently reduced. 
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1 Introduction 

It is widely believed that the early universe goes through a period of ac- 
celerated expansion called 'inflation' pQ. This believe has been successfully 
confirmed by WMAP measurements of the Cosmic Microwave Background 
radiation, which validate inflation as a mechanism to produce the initial 
irregularities of our universe [2]. Furthermore inflation has been shown to 
solve a number of cosmological problems, such as the 'horizon problem', the 
'flatness problem' and the 'monopole problem'. There are many different 
inflationary models [3]. Experimentally these can be distinguished by their 
predictions for the main observable quantities: the density perturbation at 
decoupling, the spectral index of the CMB anisotropy and the cosmic string 
contribution to the CMB anisotropy [I]. 

These observables are derived from the precise shape of the effective 
scalar potential of a particular inflation model. Only potentials that con- 
tain a nearly flat direction give rise to inflation. Supersymmetric models 
contain these nearly flat directions naturally (without fine-tuning parame- 
ters), because of large cancellations of bosonic and fermionic 1-loop quantum 
contributions to the effective potential, which give a slight tilt to classically 
flat directions. In supersymmetric hybrid inflation models, inflation occurs 
while the scalar fields are on the so called 'inflationary valley' and it ends 
abruptly with a period of spontaneous symmetry breaking. On the infla- 
tionary valley the 'Coleman- Weinberg' formula [5] can be applied to get the 
1-loop quantum correction to the effective potential [§J. This gives the shape 
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of the potential from which the aforementioned observables can be derived. 
For this derivation to be valid, it is crucial that inflation ends abruptly at the 
end of the inflationary valley, when spontaneous symmetry breaking occurs. 
This is usually assumed. 

To make precise statements about this assumption, one needs to know 
the 1-loop corrections outside the inflationary valley, during the spontaneous 
symmetry breaking phase. We will show that one cannot simply generalise 
the 'Coleman- Weinberg' formula to this case, because this introduces a non- 
renormalizable cut-off dependence and moreover it gives a 1-loop correction 
on the inflationary valley at odds with earlier calculations. The cause of 
these problems is that the 'Goldstone' boson which enters the Higgs mecha- 
nism is not massless anymore for field space points outside the inflationary 
valley. If we include this 'Goldstone' mass in the calculation of the 1-loop 
corrections, the aforementioned problems disappear. 

Using these off-valley 1-loop corrections to the D-term hybrid inflation 
potential, we show that it is in general not allowed to neglect inflation after 
spontaneous symmetry breaking. Bearing this in mind, a new analysis is 
made of all relevant observables, which are in good agreement with the 
WMAP data for a large portion of the D-term parameter space. 

This paper is organised as follows. Section [2] explains the basic ideas 
behind D-term hybrid inflation. The body of the paper is in section [3l It 
contains a way to calculate the 1-loop quantum corrections to the full D-term 
potential, which is essential to a better understanding of its implications for 
WMAP observables. Section lBTTI introduces the 'Coleman- Weinberg formula' 
for the 1-loop corrections. Then section 13.21 goes on with the calculation 
of the field dependent masses, needed for the Coleman- Weinberg formula. 
Section 13.31 explains a variety of subtleties in applying this formula and 
stumbles upon fundamental problems in calculating the 1-loop corrections. 
Section ?? hopes to resolve these problems and finally in section 13.41 the 
result for the 1-loop corrections is given. Section explores the usefulness 
of these 1-loop corrections. We will see that there is a possibility of inflation 
even after the spontaneous symmetry breaking has started. In chapter [5] 
we apply the 1-loop corrections to get new cosmological constraints on the 
D-term hybrid inflation model. If you are only interested in new parameter 
bounds for D-term inflation, consult section [5j This article will end with a 
small conclusion in section EJ 

2 D-term inflation 

The simplest possible superpotential that can generate D-term inflation is 

W = AS$ + $_ (1) 
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The super fields and <£_ have charge 1 and -1 under a U{V)fi g au ge 
symmetry. S is uncharged and A is a coupling parameter. The corresponding 
tree-level effective scalar potential is: 

2 

Vtree = >?\S\ 2 (|$+| 2 + |$-| 2 ) + A 2 |$+| 2 |$_| 2 + 9 — (|$+| 2 - |$_| 2 +£) 2 

(2) 

S is called the inflaton-field, £ is the Fayet-Illiopoulos gauge term (taken to 
be positive) and 5 is the £/(1)fz gauge coupling. <I> + = <£_ = represents 
a flat valley of local minima, because when both $ + and vanish, the 
tree-level potential does not depend on S. In this case the potential energy 
density equals: 

2 f 2 

V = 9 -j- (3) 

This is the false vacuum energy density, which drives D-term inflation. This 
inflationary solution (<3? + = <!>_ = 0) can be either stable or unstable, de- 
pending on the value of the 5-field. 

For 1 5*1 > S c the <£_ -direction corresponds to a (stable) local minimum of 
the potential, but for |5| < S c the $_-direction corresponds to an (unstable) 
local maximum, which will cause spontaneous symmetry breaking. The 
fields will settle down in the true global minimum (V = 0, S = 0, <]?+ = 0, 
|$>_| = y^) and inflation will stop. The $ + -direction cannot cause a similar 
effect, because it always corresponds to a (stable) local minimum. The 
critical value of the inflaton-field is given by: 

S c = 9 -f (4) 

At tree-level (classically) there is no term in the potential which drives 
the fields toward the true global minimum. Thus assuming that the fields 
start out with random values after the Planck-era, they will either settle 
down in the global minimum (V = 0, 5 = 0, 3>+ = 0, |<3?_| = with no 

appreciable inflation, or in the local minimum (V = 2~> 1^1 ^ *^ c ' ^+ = ^' 
=0), in which case inflation goes on forever. 

However, as shown in Ref. |6j, (l-loop)-quantum-corrections to the tree- 
level potential will add a slight slope to the D-flat direction^. Now the typical 



1 Quantum- mechanically the effective potential can be interpreted as the energy-density 
of the quantum-state which minimises this energy-density, subject to the condition that 
the field expectation values (in this case < $± >, < S >) are as given. So for every point 
in field-space you can define an effective potential value. There are however subtleties in 
the interpretation of the effective potential in case of a spontaneous symmetry breaking. 
These will be discussed in section (|B)l . 

2 Loop corrections only contribute if supersymmetry is broken. If supersymmetry is 
not broken, the bosonic and fermionic contributions to the loop correction will cancel, so 
that the total loop-correction vanishes. For this model supersymmetry is broken during 
inflation and then it is restored at the end of inflation in the global minimum. 
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evolution of the fields after the Planck-era will consist of three stages: I) First 
the and <£_ fields will quickly roll down to the D-flat valley (\S\ > S c , 
= 0, <£_ = 0). II) Then a period of slow-roll-inflation starts. The S-field 
rolls down until S c is reached. Ill) Now the spontaneous symmetry breaking 
puts a natural halt to the inflationary phase and the fields settle in the true 
global minimum, see figure 




Fi gure 1: This figure depicts the different stages for a typicai time evohition of the 
fieids in D-term inflation. In region I the $-fieIds quickly roll down to zero. Then slow- 
roll inflation starts in region II (also called the inflationary valley). Finally in region 
III spontaneous symmetry breaking occurs and the fields roll down to the true global 
minimum, whereby $_ attains a vacuum expectation value. 

For inflation-model builders, the most important part of the potential is 
the valley (\S\ > S c , <I> + = 0, = 0), which causes inflation (region II). 
Nevertheless the other parts of the potential can also be of interest. Region 
I is of interest if you want to study the initial conditions after the Planck- 
era [10\ lllj . In fact not all homogeneous initial conditions give inflation, 
but most will, because the tree-level potential is very steep in the direction 
toward the inflationary valley. Region III is of interest if you want to study 
reheating. The initial value problem will probably not lead to observable 
consequences, because everything happening more than about 60 e-folds be- 
fore the end of inflation will not have entered our observable universe today 
(it will have left its imprint on scales larger than our observable universe). 
On the other hand all processes occurring in region III will be of observa- 
tional interest, because they involve length-scales that are within our past 
light-cone. 

The study of quantum-corrections to the effective potential in these two 
regimes however is still far from finished. Actually the only part for which we 
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know the 1-loop-corrections for certain is the inflationary valley (15*1 > S c , 
= 0, <I> = 0) (See Ref. [6]). Work has been done however in the 
preheating region, Ref. [12} 113] . In the following we will try to generalise 
their approach and to point out possible errors. 

3 1-loop corrections 

3.1 Coleman Weinberg formula for calculating 1-loop cor- 
rections to the effective potential 

To compute the 1-loop correction to the effective potential people normally 
use the Coleman- Weinberg formula. It is the total contribution to the effec- 
tive potential of all Feynman diagrams, which have truncated external scalar 
propagators with zero momentum and one internal loop with any possible 
particle going around (in this case one of the scalar particles, a fermionic 
super-partner, a gaugefield or a gaugino). For the derivation in a general 
(non-SUSY) context see Ref. [5]. It is given by: 

^^E(-i)^<"(f) . 09 

with rrii being the mass of a given particle, the sum goes over all particles, F 
can be taken for bosons and 1 for fermions and A m is the renormalization 
mass. But this formula is really a simplification of the more general [14] : 

AF = eh str ( M "> A '' n {f) + 
_i_ Sjr(M2)A , + _J_ S(r ( Art „(^)) + _ (6) 

with A c being a momentum cut-off and \i the scale parameter. The dots 
stand for A c -independent contributions. Str is called the supertrace and is 
defined as follows: 

Str (M n ) = ("l) 2Jl (2Ji + 1) m? , (7) 

i 

with rrii being the field dependent mass-eigenvalues, Jj being the spin of the 
corresponding particles. The supertrace weights the number of degrees of 
freedom and gives fermions the opposite sign of bosons. 

It is important to emphasise that the two A's are not the same. They 
are sometimes confused because they appear in the same form (that's why 
there is an added subscript). The renormalization mass A m is an arbitrary 
mass-scale the theory should not dependent on, which is guaranteed via the 
renormalization group equations. The momentum cut-off A c on the other 
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hand will go to infinity after renormalization, so the potential had better 
not depend on this after renormalization. Both A's appear in a similar way 
in the formulae. In fact equation ([5]) comes from the last term in equation 
([6]) after renormalizing (which makes A c -dependence of this term disappear, 
but introduces a 'nicer' A m -dependence, which happens to be of the same 
form). 

Under what circumstances can we use equation @ rather than equation 
([6])? The first term in ([6]) will drop out in any super symmetric theory, be- 
cause Str(M°) equals the number of bosonic degrees of freedom minus the 
number of fermionic degrees of freedom, which will be zero in any super- 
symmetric theory. What about the second term? In Ref. [15] it is derived 
that under suitable requirements (meaning no anomalies) for supersymmet- 
ric theories this term also vanishes. If however we look at the derivation 
of this cancellation we see that explicit use has been made of the fact that 
all first order partial derivatives of the tree-level effective potential to the 
fields vanish. This shrinks the region of applicability of the simple Coleman- 
Weinberg formula. Only at extremum (or saddle) points are we justified to 
throw away the second term of equation ©. For the case at hand this 
means that we can only use equation ([5]) in the global minimum and along 
the inflationary valley, because here all derivatives in the tree-level potential 
vanish. 

We stress this point, because these formulae are not always used with the 
appropriate care. In Ref. [13] the 1-loop corrections to the reheating-part of 
the potential are calculated, using the simplified formula, without checking 
whether the other part really drops out. Let us now look at this case in full 
detail. 



3.2 Calculation of the field dependent masses in D-term in- 
flation 

3.2.1 Introduction 

To calculate the 1-loop correction we first have to find the masses of the 
fundamental particles in the theory. These will be dependent on the field- 
values and so in the different regimes we may have different masses. Mass- 
terms are all terms in Lagrangian that are of second order in the fields, if 
we expand around any given point in field-space. The most appropriate 
expansion is as follows: 

S = -j=(S 1 + s l + i(S 2 + s 2 )) (8) 

$ + = -^($ + l+</>+l+i($ + 2 + </>+ 2 )) (9) 

= + ^y (e+e) . (io) 

v2 
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All capital letters denote the fixed values of the fields around which we make 
an expansion. The small letters can be interpreted as the field-excitations 
with a certain mass. Of course other choices can be made, for example 
and <3?_ are treated differently. 

The well known expansion around the inflationary valley consists of set- 
ting $+1, $+2; $ and to zero and taking |5| > S c . If we want to know 
the loop-corrections off-valley in principle all fixed values are nonzero. I will 
comment later on this possibility, but for the moment it suffices to say that 
this will give very messy calculations, so we should be more restrictive in 
our choice of fixed field values. In fact we can set <£+i and <&+2 to zero, if we 
restrict our attention to regions II and III, because the spontaneous sym- 
metry breaking part for the D-term potential will be in the $_-sector. As 
discussed before the <3>+ will always be in a stable local maximum <3? + = 0. 
One could also argue that nothing should depend on the phases of the S- 
and <E>_-fields, since the tree-level potential is also not dependent on these 
phased). So one could take 5- and ^--fields that are real. Here we will take 
them complex from the start, but we will indeed see in the end that the 
1-loop correction is independent on these phases. 



3.2.2 The scalar masses 

Plugging the expansion from equations (18|9|10p into the effective potential, 
equation ([2]), we can immediately read of the mass of the 4> + \ and 4> + 2 fields. 
For the si, S2 and <ft fields we first have to diagonalize the second order part 
of the Lagrangian in order to find the physical particles (which should have 
simple propagators, so no cross-terms): 

, / ±A 2 $ 2 A 2 $5x \ / ai \ 

± ( 81 s 2 ) ±A 2 d> 2 A 2 $S 2 *2 

2 A 2 d>S 2 \\ 2 {S 2 + S 2 2 )-g 2 t + lg 2 <S> 2 ) \ tj> J 

(11) 

The corresponding (mass) 2 -eigenvalues and (non-normalised) eigenstates 
are given in table (pQ), where we used the following definitions: 

B = \\ 2 (Sl + Sl + ^ 2 ) + ^g 2 ^ 2 -g 2 i (12) 

D = \\ 2 (S 2 + S 2 -^) + \g 2 ^ 2 -g\ (13) 

C = h^ 2 (S 2 + S 2 ) + \\ 2 g 2 ^ 2 - \\ 2 g 2 & • (14) 
Now we can give a complete list of all scalar-masses in the theory, see table 

3 unlike in the F-term inflation case, where the tree-level potential does depend on the 
phases of the $+- and $_-fields. 
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(mass) 2 -eigenvalue 


(non-normalised) eigenstate 


±A 2 $ 2 


^( o ) 


i (b _ V5 2 + 4c) 


/ 2A 2 $5i \ 

2A 2 $S 2 
\D- V-B 2 + 4C / 


± (b + V^ 2 + 4C) 


/ 2A 2 $Si \ 

2A 2 $S 2 
\ D + yjB 2 + AC ) 



Table 1: Diagonalization of the (si, S2, </>)-sector 



field 


d.o.f. 


(mass)' 2 


0+1 


1 




0+2 


1 


...same... 


lin. comb, of s\, s 2 


1 


^A 2 * 2 


lin. comb. ofsi,S2, 


1 


i V£? 2 + 4CJ 


lin. comb. ofsi,.s 2 ,0 


1 


± (B + VS 2 + 4CJ 



Table 2: Mass and degrees of freedom of the scalar fields 
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3.2.3 The gauge boson mass 

In the previous section we did not get a mass for the #-field. Actually this is 
because we have chosen to write the complex field <3?_ in (jlOH in terms of a 
modulus <j) an d an angular part 6. This #-field has a nontrivial kinetic term, 
which does not allow us to interpret the second-order-part as a mass, like 
we did for the other fields. However the reason for writing things out like 
this is that the interaction of the 0-field with the massless gauge boson 
will cancel all ^-dependence of the Lagrangian in favour of a gauge boson 
B^, which acquires mass. This is the so called Higgs-mechanism. 

<I> has a charge —1 under the local gauge transformation^. Hence the 
gauge-invariant kinetic term for <I> is: 

= = (0„ - i<7A M )*l(3" + i g A»)*- 

= X -d^4> + i(* + 4>) 2 [dufldPO + g 2 A ll A» + 2gdp6A») . (15) 

There is a 'nasty' kinetic term for 9 and an unwanted coupling of d^O to A 11 . 
This quadratic coupling means we have not identified the physical particles 
correctly (it represents a mixed propagator). What we should do is write 
everything in terms of the gauge-invariant combination = A^ + ^d^6. In 
doing so, equation (fT5|) gets a simple form: 

J2W«_ = \%W4> + \g 2 {® + <t>?B^ 1 . (16) 

All ^-dependence has cancelled and the gauge field has acquired a mass (a 
mass-term for a gauge boson looks like ^rn^A^A^). The massive -B^-field 
now has three degrees of freedom instead of two for the massless A^-field. 
The #-field behaves as the longitudinal component of B^. The kinetic term 
for the gauge field remains unchanged under the redefinition, because the 
partial derivatives in F„ v = (d^A u — d u A^) arising from the shift commute. 
The mass of the gauge field is included in table © • 



field 


d.o.f. 


(mass)' 2 




3 





Table 3: Mass and degrees of freedom of the gauge field 

In principle the <£ + -field should be treated in the same way, but in this 
case it gives just simple kinetic- and mass-terms, because we expanded 
around the point <I> + = 0. If you expand around another point, you will 
get a more complicated interaction, but it will still be possible to define a 
Higgs-mechanism which works for all field-space points. 

4 under the gauge transformation $ — > e' q ^& and — > Ay, + ^<9 M £, with q the charge 
of the fields, so q = ±1, for $±,5. 
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Note that this derivation is done at any field value so not only in 
the global minimum, as is done usually. The striking thing is that the 6- 
dependence really drops out everywhere. A possible exception being the 
origin = 0, because here there is a coordinate singularity, which makes 



ill defined. 



3.2.4 The fermion masses 

In any supersymmetric theory all scalar bosons have fermionic super-partners 
and all gauge- fields have corresponding fermionic gaugino-fields, see for ex- 
ample Ref . |16j . The terms quadratic in the fermionic fields are summarised 
in the following matrix, which must again be diagonalized in order to find 
the field-dependent fermionic masses of the theory: 



( Va ip- V>+ i>s ) 



/ o 

g<f>e~ i0 


V 



A 

v/2 



g<S>e ie 


(Si - iS 2 








V2 










\ 



/ 



(17) 

ip- , ip+ and ips are the Weyl-fermions belonging to the chiral multiplet and 
tpA is the gaugino. The corresponding mass-eigenvalues are: 



1 



with 



± -\/E± ^E 2 -8X 2 g 2 <P 4 , 



E = X 2 (Sf + Sl + $ 2 ) + 2g 2 <$> 2 



(18) 



(19) 



This completes the calculation of the fermionic masses. The results are 
summarised in table 



field 


d.o.f. 


(mass) 2 


lin. comb, of ipA, ~4>-, ip+, ips 


4 


\(E + y/E 2 - 8X 2 g 2 <S> 4 j 




lin. comb, of ipA, ip-, ip+, ips 


4 


\ [E - VE 2 - 8X 2 g 2 <S>^ 





( i> A \ 
V ips J 



Table 4: Mass and degrees of freedom of the fermionic fields 



3.3 Calculating the 1-loop corrections to the D-term poten- 
tial 

Insertion of all masses in equation ([5]) in order to get a formula for the 
total 1-loop correction is in principle straightforward. Some details can be 
found in the corresponding appendices, first of all massless particles do 
not contribute (Appendix [X]). For tachyonic masses we should only take the 
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real part of the corresponding 1-loop correction (Appendix[B]) and imaginary 
mass-squared values do not occur (Appendix |C"|) . 



3.3.1 Cut-off dependence 

Since we have obtained well defined masses for all particles as a function of 
Si , S2 and <I> we can explicitly calculate the second term in equation ((6|) : 

' Str (M 2 ) A 2 C . (20) 



32tt 2 



The supertrace should be equal to zero in order to make sense of the theory. 
Otherwise, when A c goes to infinity (or to a very large value like the Planck- 
mass), the whole term blows up. Indeed in calculating the supertrace most 
contributions cancel out, but not all. We're left with: 

Str(M 2 ) = -\\ 2 (S 2 + S 2 ) - \g 2 <$> 2 + g\ (21) 

There is no renormalization-counterterm to cancel this field-dependent di- 
vergence. So we conclude that the effective potential is really ultraviolet 
divergent, except possibly for field-values where equation (|2Tj) vanishes. 
This happens exactly along the "spontaneous symmetry breaking valley", 
by which I mean the valley determined by minimising for a given Si and 
S 2 : 

= =► \\\S 2 + S 2 ) - \g 2 & + g 2 i = V * = . (22) 
3.3.2 Limit to the inflationary valley 

A good check as to whether the method in the enlarged parameter-space 
makes sense is to check if it gives the same results as before on the inflation- 
ary valley ($ = 0). Mathematically you could argue that there is a problem 
in taking the limit $ — * 0, because B„ is not well-defined for = 0. If you 
insist, take a very small value for then will be well defined and the 
Lagrangian will have no dependence on 6. The particle- masses we get along 
the valley are summarised in table ([5]). 

The 1-loop corrections on the inflationary valley are well known and easy 
to calculate if we just set <E>_ equal to zero from the very beginning. The 
particle- masses that contribute are given in table ([6]). 

Comparison of these two tables points out that the two methods really 
give different results. In the usual case there is an extra scalar particle of 
(mass) 2 equal to iA^Sf+Sf) - 5 2 £- As a result the 1-loop correction derived 
from both methods will also differ. This makes the enlarged parameter-space 
method not trustable. 
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field 


d.o.f. 


(mass) 2 


0+i 


1 




0+2 


1 


...same... 




1 







1 








1 


$X'(STi + S$) - ft 




2 or 3 





lin. comb, of ipA, V'-j VM-> V'S 


4 




lin. comb, of V'A , V'- > V'+ j V'S 


4 






Table 5: Mass and degrees of freedom of all fields along the inflationary valley in the limit 
<& — > of the enlarged parameter-space method. The gauge field has 2 degrees of freedom 
if 3> is really zero and three degrees of freedom if 3> is very small. For the computation of 
the 1-loop corrections it makes no difference (2-0 = 3-0 = 0). 



field 


d.o.f. 


(mass)' 2 


0+1 


1 


LX\Sl + Sl)+g^ 


0+2 


1 


...same... 


0-1 


1 


+ 


0-2 


1 


...same... 


lin. comb. ofip—,i/)+ 


4 





Table 6: Mass and degrees of freedom of all massive fields along the inflationary valley 
as it is commonly used. 
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3.3.3 The cut-off dependence and the limit problem 

Obviously the described method similar to the one in Ref. [13] cannot 
be correct. Results are unphysical because of the cut-off dependence and 
just wrong because of the contradiction with the previous method along the 
inflationary valley. In the following we will describe a possible method which 
doesn't suffer from these two problems. 

The origin of the limit-problem seems to be the Higgs-mechanism. The 
Higgs- mechanism dictates that we write the $_ -field in polar coordinates. 
In doing so the angular #-field becomes massless. This could well be the 
mass we're missing with respect to the usual case. Furthermore the cut-off 
dependent supertrace contribution in equation (f2lT) also amounts to the same 
(mass) 2 -value of ^A 2 (5 2 + 5|) — g 2 £ in the limit <I> goes to zero. So maybe 
we can solve both problems at a time by evaluating the loop-contributions 
with the masses given by the theory before applying the Higgs-mechanism. 
First I will show that this indeed solves the aforementioned problems. Then 
we will argue about the physical validity of this method. 



3.3.4 A possible solution 

Since we suspect that the problems come from the massless #-field we will 
use the following expansion: 

S = -L(S l + s l + i(S 2 + s 2 )) (23) 

*+ = ^(^+1 + ^+2) (24) 

$- = -^(*-i + ^-i + i($- 2 + 0-2)) (25) 

Also we will set S 2 and $_2 to zero, because otherwise the calculations 
become unnecessary bothersome and we have seen before that these phases 
really make no difference in D-term inflation. Going through the same steps 
as before in calculating the scalar masses (see section 13.2.21) we get the 
following mixing between the (si, s 2 , <p-i, (fi- 2 )-Helds: 

^ ( 0_! 0-2 si s 2 ) (26) 

/ 1^+3^2^ _^ A 2 5l$_l \ 

\\ 2 Sl + \g 2 ^\-g\ 

A 2 Si$_i " \>?®-\ 

\ \\ 2 ^\ j 

Diagonalizing this (mass 2 )-matrix and going back to the original notation 
(replacing S\ by (S\ + 52) and 3>_i by <£) gives the following field dependent 
particle masses, see table (|7|). 



/ 4>-i 

<P-2 
Si 

V s 2 
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field (for S 2 = and = 0) 


d.o.f. 


(mass) 2 


<t>+l 


1 




<P+2 


1 


...same... 


S2 


1 


±A 2 $ 2 


4>-2 


1 




lin. comb. ofsi,0_i 


1 


1(b-VB 2 + 4C) 


lin. comb, ofsi,0_i 


1 


i (fi + VB 2 + 4CJ 



Table 7: Mass and degrees of freedom of the scalar fields, without applying the Higgs- 
mechanism (meaning there are still quadratic mixing-terms with left). The fields 
are given for 52 and $_2 set to zero. For other fixed values the (mass) 2 -eigenfields are 
different, but the mass-spectrum is the same. The definition of B and C is as before, see 
equations (O and (fTl)) . 



The mass-spectrum is the same as in table ([2]) with an extra mass for 
the (/>_2-field. Comparing this extra mass with equation (|21|) reveals that 
using this method the A 2 -dependence really drops out. Also In the limit 
<3? —* the scalar masses in tables (|7|) and © become identical, resulting in 
a correct smooth limit to the known 1-loop corrections along the inflationary 
valley. To be precise, in table ([7]) </>_2 and one of the linear combinations of 
(si, <p-i) get the same mass, the other linear combination becomes massless. 

This new method works because the masses of the gauge field and the 
fermionic fields remain unaffected by this different choice of coordinates, so 
all other masses remain the same as before. For the fermionic part this 
is trivial, because the fermionic masses only depend on the fixed values of 
the scalar fields. For the gauge field this is less trivial, but comparison of 
equations (fT5|) and (fT6l) reveals that the mass of the gauge field is unaffected 
by the Higgs mechanism. In both cases the (mass) 2 -value equals g 2 <3? 2 . 
In principle the difference between both methods lies in the inclusion of 
Feynman-graphs with the mixed 'Goldstone'-gauge propagator. In the next 
section we will present another way of viewing this situation. 

3.3.5 Validity of the method 

Of course we can't just change the calculational method, because there is 
a physical background behind this. The reason for trying this out anyway 
is that there is no clearcut answer in the literature which relates the Higgs- 
mechanism to the Coleman- Weinberg formula and which states explicitly 
what are the masses to be used in applying the Coleman- Weinberg formula 
in this setting. This is because normally if you only look at the global 
minimum it really makes no difference, since here the Goldstone mode is 
massless anyway and, as is shown in appendix [XJ massless fields do not 
contribute to the 1-loop correction. 
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The Coleman- Weinberg formula, as used before, includes all 1-loop Feyn- 
man graph contributions to the effective potential, where the external trun- 
cated scalar legs have zero momentum (since the vacuum has no momen- 
tum). The particle in the internal loop can be any particle in the theory 
that couples to the scalar external legs. In the cut-off dependent Str(M 2 )- 
part of the Coleman- Weinberg formula it really makes no difference whether 
we take the (mass) 2 -values of the physical particles, after diagonalizing the 
(mass) 2 -matrix, or we take the 'unphysical' particle masses (without diago- 
nalizing this matrix) . This is because the trace of a matrix is invariant under 
a change of basis lj. Somehow in mixing the scalar field with the gauge field 
via the Higgs mechanism, it seems that it does make a difference which basis 
we use. 

In appendix |D] we describe a Higgs mechanism for a simple U(l) invari- 
ant potential, using both polar and Cartesian field definitions. Surprisingly 
for the global U(l) case the 'Goldstone' boson has a mass, if it is outside the 
global minimum. For the local case one usually uses polar fields, because 
of the Higgs mechanism. However it is also possible to define the equiv- 
alent Higgs mechanism with Cartesian fields. Then the mixed propagator 
term vanishes, as does the kinetic term for the 'Goldstone' field, but the 
'Goldstone' mass-term does not vanish^. 

Of course the discussion in appendix|D]holds equally well for the Mexican 
hat type U(l) symmetry we're interested in. In this case the (mass) 2 -value 
of the global Goldstone mode will be a smooth function of <£, which goes from 
a negative value at the top of the Mexican hat smoothly to zero at the global 
minimum circle to positive values outside this circle, see figure ([2|). Surpris- 
ingly the (mass) 2 of this Goldstone mode ^X 2 (S 2 + Sf) + \g 2 ^ 2 — g 2 i is 
exactly what we need to solve the cut-off dependence and the limit problem. 
So using second derivatives of the potential for the classical masses in the 
Coleman Weinberg formula gives a well-defined 1-loop correction even out- 
side the global minimum and using the tree-level propagator masses doesn't. 

There are now two ways of viewing the situation. 

1. We should not include the Goldstone masses in a calculation of the 
1-loop corrections to the effective D-term potential outside the global 
minimum. As a result these 1-loop corrections are ill-defined for almost 
all field values. 

2. We should include the Goldstone masses in the calculation of the 1- 
loop corrections to the effective D-term potential outside the global 

5 Use Tr(AB) = Tr(BA). If matrix M is diagonalized as follows M = CDC' 1 , then 
Tr(M) = TriCDC' 1 ) = Tr^CC^D) = Tr(D). Also Tr(M n ) = Tr(D n ). 

6 Maybe 'mass-term' is not the right word to use here, because for the local U(l) case 
it no longer corresponds to a pole in the propagator. It is the term quadratic in the 
corresponding 'Goldstone' field 
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Figure 2: The solid line gives the (Mass) 2 -value of the global 'Goldstone' mode as a 
smooth function of the <3?-field for S = 0. The dashed line depicts the shape of the tree- 
level effective potential (the vertical axis doesn't refer to this). Note that the Goldstone 
mode is indeed massless at the global minimum of the potential. 



minimum. The resulting 1-loop corrections agree with previous results 
along the inflationary valley and have no A^-dependence. This is the 
viewpoint we will take. 

There are some field-space points for which both methods agree, because 
the Goldstone boson is really massless there. These points lie exactly on the 
"spontaneous symmetry breaking valley" defined by equation (I22D . Note 
that both methods disagree on the inflationary valley. 



3.4 Graphs of the one-loop correction 

This section is merely meant to give the reader a feeling of what the 1-loop 
corrections look like in comparison with the tree-level potential and what 
are the possible implications. A quantitative analysis of the different aspects 
will be given later on. The results are shown with the aid of graphs, because 
the actual algebraic expression is too complicated^- 

There are four adjustable parameters: A, g, £ and A m , although the 
dependence on the last will be small. The output, being a potential energy 
density, will have mass dimension four and will always be given in units of 
(Mpi) A . Don't get confused by very small numbers, this is just because of 

7 In fact the slightly broken supersymmetric nature of the 1-loop correction makes 
it of the form: near-infinity minus near-infinity. Huge bosonic and fermionic contribu- 
tions nearly cancel and the difference is the correction we are interested in. Many of 
these cancellations can be done analytically, by tricks of the form: A 2 ln(A) — B 2 ln(B) — 
A 2 ln(A/B) + (A 2 - B 2 )ln(B), which can make things much simpler if A and B are, 
nearly equal, huge numbers or formulae. Unfortunately we didn't manage to get rid of 
all (near-infinity minus near-infinity)-terms in this way, so the final formula for the total 
1-loop correction is still too complicated to analyse by hand. The near-infinity minus 
near-infinity nature can also lead to numerical problems, but these we have solved. 
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the fourth power. First we will look at the case where all these parameters 
are one, although this will not be the typical case, because normally £ is 
much less than one square Planck mass. 

Figure ©(a) gives the classical tree-level potential in the region < S < 
2 • S c . The <3?-axis reaches to 1.2 times the global minimum value, so we can 
see the entire classical path during spontaneous symmetry breaking. Part 
(b) depicts the 1-loop correction on the same domain. It is around a factor 
twenty smaller than the tree-level potential and has an opposite shape. Of 
course we rely on the 1-loop corrections to be (very) small in comparison 
to the tree-level potential, so that we are hopefully justified to ignore all 
two and higher order corrections. Part (c) depicts the total (tree + 1-loop) 
potential. As mentioned in section ([2]) the quantum-correction indeed adds 
a little slope to the inflationary valley ($ = 0). In slowly rolling down from 
2 • S c to S c , this universe will undergo around 2900 efolds. 

The 1-loop corrections in this case are small as compared to the tree level 
potential. In most cases they will be even much smaller, but this doesn't 
mean that they are unimportant, because sometimes it is the variation that 
matters. In figure ([!]), we have zoomed in on the tree-, 1-loop- and total 
values of the potential in a region centred around the critical point (S c ). This 
part of the inflationary valley corresponds to around 52 efolds. Although the 
tree-level potential is a factor 60 larger, the 1-loop corrections still dominate 
the behaviour around the critical point. You can nicely see the slope that 
the 1-loop corrections generate. The assumption that inflation ends when 
the fields leave the inflationary valley, which is usually made, is probably 
not valid here: since the upper half of the picture corresponds to 52 efolds, 
we can already see that the total potential doesn't change fast enough on 
the relevant time-scales to end inflation abruptly (note that the scales on 
both axes are the same). The dependence on the renormalization mass scale 
(A m ) is very small. If we change the value of A m from 1 Mpi to 0.01 Mpi, 
the slope changes only to give around 58 efolds instead of 52 on the same 
domain. 

Figures © and (jBJ) represent the general shape of the 1-loop corrections, 
which you get for almost all physically relevant parameter combinations. 
For these figures g = 0.1, A = 0.0001, ^ = 0.00025 M Pl and A m = 0.35 
Mpi (which is equal to the value of S c ). These values seem a bit out of 
the blue, but in fact I took a combination which is compatible with WMAP 
observations, which will be explained in section ©. Anyway, this selection 
doesn't change the overall picture. 

The main difference with the previous set of figures is that the 1-loop 
corrections are even smaller as compared to the tree-level potential (around 
a factor of 500) and that the 1-loop corrections are only of opposite shape 
as compared to the tree level potential for small $>. For large $ the 1-loop 
corrections mimic the tree-level potential. So only in the neighbourhood of 
the inflationary valley can we expect the corrections to be of interest. 



18 




(a) tree-level (b) 1-loop (c) total 

Figure 3: The tree level potential, 1-loop correction and the sum of both, for g — 1.0, 
A = 1.0, vf = 1.0 Mpi, A m = 1.0 Mpi on the domain < S < 2- S c , < $ < 1.2-$ g . min . 
The tree-level potential looks similar to figure {1} . The 1-loop correction is around a factor 
20 smaller than the tree-level potential and it has an opposite shape. The general shape 
is not changed by the 1-loop corrections, but there is a slight slope on the inflationary 
valley, which can not really be seen on this scale. It is very small, which results in many 
efolds: around 2900 in the upper half of the picture. 




(a) tree-level (b) 1-loop (c) total 

Figure 4: The tree level potential, 1-loop correction and the sum of both, for the same 
parameters as in figure © on a domain centred around the critical point (S c ). Both field- 
axes have the same scale. The slope along the inflationary valley, arising from the 1-loop 
corrections is clearly visible now. It will give rise to 52 efolds of inflation on the upper 
half of the picture, which therefore represents the cosmologically interesting region. It is 
apparent that the variations in the 1-loop corrections dominate the tree-level variations 
in this region. 



Another thing, which is shared by all possible parameter choices, is that 
in the global minimum both the tree-level potential and the 1-loop correc- 
tions vanish, giving a true, totally supersymmetric, non inflating, ground 
state. 

In the (S = 0)-plane the 1-loop corrections always counteract the tree- 
level potential, although normally the effect is too small to have important 
consequences. 

Close to the inflationary valley the 1-loop corrections behave more like, 
say, a parabola and the tree-level potential more like a fourth order function. 
So the closer you get to the inflationary valley the more important relatively 
the 1-loop corrections become. 

Again, if we zoom in on the relevant region around the critical point 
we see that the 1-loop corrections here will definitely have important con- 
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sequences. This opens an interesting window. It seems that the, up to now 
neglected, off-valley quantum-corrections will play a major role on all scales 
which are of cosmological interest. Most surely they will have consequences 
to the explanation of WMAF@ observations. 




(a) tree-level (b) 1-loop (c) total 

Figure 5: The tree level potential, 1-loop correction and the sum of both, for g = 0.1, 
A = 1.0 • 10~ 4 , = 2.46 ■ 10~ 4 Mpi and A m = 0.35 Mpi on the domain < S < 2 ■ S c , 
< $ < 1.2 ■ & g .min.- The 1-loop correction is around a factor of 500 smaller than the 
tree-level potential. The shape is opposite near the inflationary valley. The tilt produced 
by the 1-loop corrections corresponds to around 18000 efolds in the upper half of the 
picture. 




(a) tree-level (b) 1-loop (c) total 

Figure 6: The tree level potential, 1-loop correction and the sum of both, for the same 
parameters as in figure (JSJ zoomed in around S c - Both axes have a different scaling. The 
upper half of the picture corresponds to 58 efolds, so the 1-loop corrections are definitely 
important on the relevant scales. 

For a select region in parameter space it is possible that the negative 
1-loop corrections at the (S, <£) = (0,0) point become relatively large, say 
order 10% of the tree-level potential. This less common situation is depicted 
in figures ([7]) and ([8]). The lowering of the top of the Mexican hat would 
probably lead to less energetic cosmic strings. Of course if the 1-loop correc- 
tions become too big, we can no longer neglect the second and higher order 
loop contributions. Again we show a parameter combination which will be 
compatible with the WMAP density perturbation. 

8 The Wilkinson Microwave Anisotropy Probe is a satellite which looks in great de- 
tail at the cosmic microwave background radiation, see |http:/ /map.gsfc.nasa.govj It's 
predecessor was the COBE satellite. 
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(a) tree-level (b) 1-loop (c) total 

Figure 7: The tree level potential, 1-loop correction and the sum of both, for g = 1.0, 
A = 0.01, = 6.8 ■ 10~ 4 M m and A m = 0.096 M P i (= S c ) on the domain < S < S c , 
< $ < 1.2 ■ <& a .min. • The 1-loop corrections have lowered the top of the Mexican hat 
(origin in this picture) considerably and will thus slightly lower the resulting cosmic string 
tension. 




(a) tree-level (b) 1-loop (c) total 

Figure 8: The tree level potential, 1-loop correction and the sum of both, for the same 
parameters as in figure on the domain < S < 2 • S c , small $. This shows how 
the height of the Mexican hat goes down as the inflaton field S goes to zero. The 1-loop 
corrections will slightly decrease the tension of the resulting cosmic strings. 



4 Inflation below the critical inflaton value 

Usually it is assumed that no significant inflation occurs after the inflaton 
field S reaches the critical value S c . In this section we will show that even 
classically (without 1-loop corrections) this assumption may be false. Then 
we will calculate explicitly the number of e-folds below S c for the D-term 
parameter combinations which best fit WMAP measurements, according to 
Ref. |20j. We will see that indeed inflation below S c can not be neglected. 
A whole new analysis is needed in order to compare the predictions of the 
D-term model with WMAP measurements. This is done in section ([5]). 

Substantial inflation after S c is interesting for two reasons. First of all 
this offers the possibility to dilute the cosmic strings which are produced 
in the symmetry breaking process. This is a good thing, since WMAP has 
put stringent bounds on the string contribution to the cosmic background 
anisotropy. It can not be more than about 10% H pH ITS] . 

9 The cosmic microwave background anisotropies have a fluctuating nature if they arise 
from quantum fluctuations that are enlarged by inflation. This is because of acoustic 
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Secondly, the cosmologically interesting region, which is around 50 to 60 
efolds before the end of inflation, will shift when substantial inflation occurs 
after S c . This region, for example, determines the density perturbation at 
reionization (which should be around 2 • 10 -5 in order to be compatible 
with WMAP measurements). So the density perturbation at reionization is 
dependent on the possible inflation below S c . 

Inflation takes place whenever the potential energy-density of the scalar 
fields dominates the kinetic energy density [3]. This is possible only if the 
slope of the potential is very small. In this case the equation of motion 
becomes dominated by the friction term (Hdt4>), the fields are in slow-roll 
and the number of efolds before the end of inflation is effectively determined 
by how long the fields manage to stay in slow-roll. A good approximation 
to the number of efolds is [3] : 




The integral is over the path in field-space, units are in Planck masses and 
the approximation is valid as long as the slow-roll conditions are satisfied. 

One normally assumes that inflation ends either before reaching S c , 
whenever the slow-roll parameters become too big (order one), or at S c , 
because then spontaneous symmetry breaking starts and will speed up the 
fields very quickly. In the previous chapter we have seen that this is prob- 
ably not always true, because the potential does not always become very 
steep in a region around S c corresponding to about 50 to 60 efolds. Partly 
this is because close to the inflationary valley the 1-loop corrections are just 
opposite to the classical potential. But also even classically this assumption 
may be false. 

Figure (J9j) gives the general shape of the tree-level potential in terms 
of three parameters a, b, c which correspond respectively to the width and 
height of the Mexican hat and to the value of S c . The relation with the 
D-term parameters is: 

a = VM Z = \a 2 

b=^f or g = ^ . (28) 

c _ gV^ \ _ 2V2E 

A ac 

The conditions that favour inflation during spontaneous symmetry breaking 
are a large a, small b and large c, because these make the potential flat. Since 
the dependence on the physical parameters is fairly simple this translates to 
a large £, small g and small A. We can even get a simple bound on when 
inflation happens approximately all the way down to the global minimum. 

oscillations in the period between inflation and reionization. Anisotropics coming from 
cosmic strings have no specific scale to it. WMAP measurement look very much like the 
"quantum fluctuation" anisotropies, leaving at most a 10% contribution to cosmic strings. 
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Figure 9: This picture depicts the typical shape of the spontaneous symmetry breaking 
part of the classical D-term potential. Inflation during spontaneous symmetry breaking is 
favourable if the potential is flat, meaning large a, small b and large c. 



Let's assume that a is much smaller than C0 and the field follows a straight 
line in going from S c to the global minimum, then we have the following 
average values: 

<y>=-j r (29) 

< V >= 9 | ~ l±- (30) 
2Va 2 + c 2 2c V ; 

1 (<v>\ 2 i , 01 v 
<e>= i^U~^J (3) 

with e being one of the slow-roll parameters [3]. The condition for inflation 
to persist approximately all the way down becomes very simple: 

e < 1 , S c > —= « 0.28M PJ . (32) 

V47T 



Inflation will happen during a large part of the spontaneous symmetry break- 
ing if S c is bigger than 0.28 Mpi. Otherwise it is still possible that the area 
around S c is flat enough to give substantial inflation in the beginning of the 
symmetry breaking, but certainly not during this whole stage. 

We will now look at the parameter space which give the right magni- 
tude of density perturbation according to Ref. [20]. Figure 2 in Ref. |20j 



10 This is quite natural to assume, because big values of £ will give cosmic strings with 
a too large tension |18| . 
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gives various parameter combinations that lead to the right density pertur- 
bation^- 

Moreover combinations that lie under the horizontal line give a 
cosmic string contribution less than 10%. 

Using their parameters as input, we are able to see if there is any in- 
flation after S c (which is not included in their computation). We calculate 
numerically the steepest path down, using of course the total (tree-level + 
1-loop) potential and, with the aid of equation (|27|) . we attain the number 
of efolds after S c . Figure (fTUj) gives the results for several datapoints with 
g = 0.01. Inflation stops when either e or rj comes close to one. You can see 
that it doesn't really matter which criterion we use. In general r/ will become 
large just before e does so and even if there is a substantial difference in the 
field values at which respectively the r\ and e criterions fail, this will not give 
a substantial difference in efolds, because when rj becomes large the slope 
is already to steep to have a considerable contribution from equation (|27j) . 
However e is a more reliable numerical quantity than rj, since it is dependent 
on the first derivative of V instead of the second derivative. It is important 
to note that in all simulations the difference between both criteria never 
exceeds about 3 efolds. Since the r\ criterion can sometimes be off due to an 
enhanced numerical error, we will stick to the e criterion from now on. 
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Figure 10: This figure depicts the number of efolds after S c using datapoints from Ref. 
[20] for g — 0.01. In one simulation inflation stops when e becomes one, in the other case 
when r\ becomes one. This gives approximately the same results. 

Figure (jlip gives the results for all sampled parameter combinations. 
The remarkable thing is that the number of e-folds after S c can be as high 

11 Note that k in Ref. |20) is A in my notation. 
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Figure 11: This figure depicts the number of efolds after S c using datapoints from Ref. 
|2U] spanning the whole range of parameters. As you can see in approximately half of the 
cases the number of efolds exceeds 10, which will certainly give detectable consequences. 
In some cases the number of efolds even exceeds 60, meaning that the whole region of 
inflation relevant to cosmology happens after Se- 



as 100000, so we are definitely not always justified to ignore this. It is hard 
to extract further physical results from this, because as soon as the number 
of efolds after S c become substantial, it could well be that this messes up 
the WMAP agreement in the density perturbation, which is dependent on 
the region about 60 efolds before the end of inflation. In other words, the 
relation between A, g and £ used in Ref. [20J could be wrong. In the next 
section I will try to find the corrected relation between these parameters, 
that is compatible with WMAP. 
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5 New bounds on SUSY D-term inflation param- 
eters 



We have found a way to compute 1-loop quantum corrections to the D- 
term effective potential in section ©. In section (??) we have seen that 
these 1-loop corrections might be more important than thought previously. 
Especially in section we saw that, using these 1-loop corrections and 
looking at data from [20J , there can be a huge amount of inflation after the 
critical S- value (where inflation is assumed to stop). The purpose of this 
section is to re-examine the agreement with the WMAP data [2], given the 
new insight of inflation after spontaneous symmetry breaking. 

5.1 Matching with WMAP-data 

There are three free parameters g, A and £. g is the dimensionless gauge 
coupling. Its most natural value would be somewhat smaller than unity. A is 
a dimensionless coupling and £ is the Fayet-Illiopoulos gauge term in units of 
(Mpi) 2 . Its squareroot gives the field-value of the field after spontaneous 
symmetry breaking, see equation (|28|) , which directly determines the tension 
of the resulting cosmic strings, see equation ([39]) . The energy scale at which 
spontaneous symmetry breaking occurs is equal to \/£ up to a factor g/X, 
see equation (jlj). 

The formulae for the amplitude of the (scalar) density perturbations A$ 
and gravitational waves (or tensor perturbations) At are[?]: 



The input is again in Planck masses and both should be evaluated at 60 
efolds before the end of inflation (or later if you want to know the perturba- 
tions on smaller scales). The value of the potential energy density V and its 
derivative in the direction of steepest descent V' are coming from the total 
(tree level plus 1-loop) potential. As long as At is small in comparison with 
As it can be neglected and A$ should be approximately 2-10 -5 according to 
CMB measurements. In all simulations, the amplitude of the gravitational 
waves is very small and thus can be neglected. 

We proceed as follows: we fix a given combination of A and g, set A m 
equal to S c so that we can compare results with Ref. then we guess 

12 Normally the 1-loop corrections are very small in comparison to the tree-level poten- 
tial. That is why people don't include the 1-loop corrections in the value of V plugged in 
to equations like (I33p and (|34p . The nice thing is that derivatives of V are independent 
of A m . Neglecting the 1-loop corrections in V is practically the same as setting A m equal 




(33) 




(34) 
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a value of £. We do the whole simulation with this value of £, see whether 
we get the wanted value of 2 • 10 -5 for the density perturbation, equation 
(|33f) . If not, an algorithm makes a good second guess, and so on, until the 
density perturbation converges to the wanted limitP^l. 

The main differences between our analysis and the analysis in |20| are: 

1. Their simulation stops at S c , because they only have the 1-loop correc- 
tions on the inflationary valley. This new simulation takes into account 
the number of efolds after S c and the off- valley 1-loop corrections. 

2. Because of aforementioned numerical problems of the kind " near-infinity 
minus near- infinity" , see footnote (|7|) , in Ref . [20] they use an approx- 
imate formula for the 1-loop corrections. We solved these numerical 
difficulties and use the real formula. 

3. The simulation in Ref. [20] includes the cosmic string contribution in 
the determination of £. Our simulation doesn't. This is justified since 
we can not assume right away that cosmic strings contribute, due to 
inflation after spontaneous symmetry breaking, see section (I5.4p . If 
the cosmic strings do contribute this will give a difference of at most 
10% for datapoints consistent with WMAP. 

This now happens to be one of the rare occasions in which increasing 
the complexity of the computation actually makes the results more simple, 
as can be seen by comparing figure (fT2j) with figure 2 in Ref. [20] . It depicts 
the whole range of values of the three parameters A, g and £ that give the 
right density perturbation of 2 • 10 -5 . The parameter-range is larger than 
in Ref. [20]; There is an extra dataset for g = 0.3 and g = 1.0. Although 
a lot of computerwork is put in to get all these datapoints, the result is 
astonishingly simple. In Ref. [20J there is a large dependence on the gauge 
parameter g. This dependence has totally disappeared. Only for g = 1.0 
can we detect a slightly different path. For A-values up to about 5 • 10~ 4 the 
remaining (A,£)-dependence is perfectly described by: 

= 6.42 • 10 lfl A 1 / 3 (GeV) . (35) 

We cannot explain why a very involved simulation puts out such a nice 
precise number as 1/3. 

The value of the gauge coupling g has no influence on relation (f35|) , but 
it does determine at which point the number of off-valley efolds become im- 
portant. This can be seen in figure (|13|) . For each value of g the simulation 

to the energy-scale of interest, since this makes the 1-loop corrections very small in the 
direct neighbourhood of this energy-scale. As we saw before, the A m -dependence is very 
small anyway, so different choices of A m don't give dramatically different results. 

13 There is really no hope to do such a thing analytically, because the 1-loop corrections 
are way too complicated, and also the trajectory is not a straight line, as it is if you stick 
to values of S bigger than 5 C . 
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stops when the number of off-valley efolds comes close to 60. The value of A 
at which this occurs becomes smaller if g becomes smaller. The simulation 
stops at this point, because when the 60th efold is off-valley the determi- 
nation of the density perturbation becomes numerically very difficult and 
more importantly, my algorithm to find the desired £-value consistent with 
WMAP does not converge any more. We do not exclude the possibility that 
there are still parameter-combinations with smaller A-values compatible to 
with WMAP, but they don't lie in the neighbourhood of the other points. 
In mentioning the number of off-valley efolds I really mean the number of 
efolds after the <I>-field leaves the valley (<£ = 0). In general the 1-loop cor- 
rections first extend the steepest descent path a bit along the S'-axis and 
then spontaneous symmetry breaking start quite abruptly, see for example 
figure ©(c), but note that the scaling is different on both axes. This is the 
starting point for what I call 'off- valley inflation'. 

The reason why the g = 1 plot does not perfectly overlap with the other 
g- values is probably that for g = 1.0 already off- valley inflation becomes 
important around A = 5 • 10~ 4 , before the curve bends around the corner. 

The resulting parameter-regime that is compatible with WMAP is thus 
much larger than previously estimated. For small A values all datapoints 
remain at low £ values, which puts them safely below the 10% cosmic string 
contribution line at y 7 ? = 2.3 • 10 15 GeV. But even more important, before- 
hand the maximum value of g which still had some datapoints below the 
10% limit was g = 0.03 and this has gone up to about g = 0.3. In principle 
gauge couplings close to unity are more natural and they could arise from 
string theory. Even for gauge-couplings g bigger than 0.3, there might be 
combinations compatible with WMAP, because for all g's it is possible to get 
substantial off-valley inflation which might dilute the cosmic strings enough 
to suppress its contribution to the density perturbation at reionization (an 
effect which is not included in the 10% bound in figure (|12p ). If, for exam- 
ple, as a first guess we assume that there is one cosmic string produced on 
average per Hubble volume at the moment of symmetry breaking, then 60 
efolds of off-valley inflation would mean that there is on average about one 
cosmic string left in the whole of our observable universe. The contribution 
to the density perturbation would then be way less than 10%. 
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Figure 12: This graph shows the parameter-combinations that are compatible with the 
WMAP density perturbation. Points below the 10% cosmic string contribution line are 
also compatible with the WMAP measurements of the spectrum of the CMB-anisotropies. 
There seems to be almost no dependence on the gauge parameter g. Datapoints for 
different g overlap perfectly. The formula: = 6.42 • 10 16 A 1/3 (GeV) gives a perfect fit 
for A < 5 ■ 1(T 4 
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(c) g = 0.1 



(d) g = 0.01 • V10 
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(e) g = 0.01 
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Figure 13: The same dataset as in figure (|12p for different values of the gauge coupling 
g separately. The cut-off at the lower A-value corresponds to about 60 off-valley efolds. 
Taking smaller g-values shifts this cut-off to consecutively smaller A-values. 
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5.2 Number of e-folds 



The reason why the results in section (|5.ip have simplified so dramatically, 
is because of the inclusion of a descent calculation of the amount of infla- 
tion below the critical inflaton field value. This shifts the cosmologically 
interesting region, which is 60 efolds before the end of inflation and there- 
fore changes all observables, like the density perturbation and the spectral 
index, see section (]5.3p . The effect will prove to be the largest for big values 
of g and small values of A. let us now look at the number of e- foldings in 
more detail. 

As mentioned before in section @ the 1-loop corrections will lengthen 
the path travelled on the inflationary valley. There are now two variables 
to look at: the number of efolds after S c , something which was neglected 
before, but more interesting even is the number of efolds after the symmetry 
breaking, because this opens the possibility to dilute the cosmic strings after 
they are produced, which will be discussed in section (|5.4p . 

In order to compute the number of efolds we use the integral in equation 
(|27|) . Figures (|14|) to (|16j) show the number of e-folds as a function of g and 
A and the fitted £. Again the results are very simple. The following formulae 
give a good fit, except for g = 1.0. 

N Sc = 7.4 • 10-y- 94 A- 140 (36) 

JY twJ = 7.4.10-V- 91 A- 1 - 40 (37) 

Nsc and N va i give the number of efolds respectively after S c and after leaving 
the ($ = 0) valley. It is quite surprising that the fit works well in both the 
g and A direction and moreover both functions are very similar, with about 
a factor of 10 difference. 

In section (jlj) we argued that, even classically, inflation can be extended 
on a large portion of the trajectory during spontaneous symmetry breaking, 
if S c becomes bigger than about 0.28 Mpi, see equation (|32l) . Figure (fT7|) 
shows that this is indeed the case. Right around S c ~ 0.28Mpi, the por- 
tion of the off-valley trajectory on which inflation still takes place, becomes 
substantial. This does not depend much on any of the other parameters. 

5.3 Spectral index 

The spectral index denotes how the amplitude of the tiny fluctuations de- 
pend on the scale. A spectral index of n = 1 would correspond to a scale 
invariant perturbation. WMAP measurements reveal that this spectral in- 
dex is very close to one. The 1-year results gave a value of n = 0.99 ± 0.04. 
The accuracy improved with the 3-year results (which include CMB po- 
larisation measurements), giving a value of n = 0.961 db 0.017 [2]. So the 
spectral index is very close to one and probably a bit smaller (n = 1 is at a 
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lire 14: This graph shows the amount of inflation both after S c and after leaving the 
= 0)-valley for all datapoints. 
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lire 15: This graph shows the amount of inflation after S c for all datapoints. Ns c 
10 _3 <7 1 ' 94 A -1 ' 40 is a good fit to all datapoints except g = 1.0. 
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Figure 16: This graph shows the amount of inflation after leaving the (<f> = 0)-valley for 
all datapoints. Nsc = 7.4 • 10 _4 <? 1,91 A _1 ' 40 is a good fit to all datapoints except g — 1.0. 
The occasional point that is off the line is probably due to a slight mistake in the automatic 
recognition of the point where $ really leaves the neighbourhood of (4> = 0), which can 
be tricky. 
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Figure 17: The SySc field value at the end of inflation is plotted as a function of S c - At 
around the predicted value of 0.28 Mpi this line goes steeply down, but also on the left of 
the picture there seems to be a possibility of a very flat potential up to around half-way 
down the S-axis. Most datapoints however lie very crowded on one curve. 
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2<T deviation for the three-year results) . Bayesian statistics give a somewhat 
larger interval as compared to the error estimate above |21j . 

During slow-roll there is a simple equation for the spectral index PJ: 

n = 1 - 6e + 2rj . (38) 

In our computer-program this is evaluated between 58 and 60 efolds before 
the end of inflation. 

The results emerging from our simulations are very simple. Again there 
is almost no dependence on the gauge-parameter g. The spectral index 
equals 1.0 for all A smaller than 10~ 4 and then goes to a value of n = 0.983 
and remains constant for A > 0.01, see figure (I18D . compare with Ref. [20] . 
figure 2. The slight deviation for g = 1, from all other g- values, is exactly 
in the same A-region that was also slightly off in figure (]12p . 

The nice thing is that the model gives very precise predictions, so it can 
be falsified. In [20] for example the spectral index is equal to 1.0 on most of 
the domain and then suddenly goes down very fast, meaning that you can 
get any small value of n, but you have to fine-tune the parameters, which 
is not desirable. Our value of n = 0.983 is at a 1.3<7 deviation from WMAP 
measurements and n = 1.0 is at 2.2a deviation. So the values of A that 
are consistent with WMAP, are approximately A > 0.001, although smaller 
values are not excluded. 

The region A > 10~ 3 , which gives the best spectral index according 
to the 3-year WMAP data, does not overlap with the region consistent 
with the 10% cosmic-string contribution bound, which is approximately A < 
0.5 • 10 -4 or A < 10~ 3 if cosmic strings can be effectively diluted, see section 
(|5.4p . Of course, this can still be due to the statistical error in the WMAP- 
measurement (Bayesian statistics for WMAP do not exclude a spectral index 
equal to one) . It might also be possible to get a spectral index smaller than 
one, by adding fluctuations to the path. Now we took the steepest path 
down, which is a sensible thing to do because of the large friction term in 
the equations of motion, but since the path on the inflationary valley is not 
necessarily long, it could be that small deviations from the steepest path 
have not damped out and these could give a spectral index other than one 
(since the spectral index has a large dependence on the second derivative 
of the path). Another possible, but maybe far-fetched solution, will be 
presented in section (|5.5j) . 

5.4 Diluting cosmic strings 

In section (|5.2p we have seen that a large amount of off-valley inflation is 
possible for any value of the gauge coupling g. Cosmic strings will form at the 
end of inflation, because of the spontaneous breaking of a U(l) symmetry 
[9]. In this section we will see how off- valley inflation affects the cosmic 
string density. 
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Figure 18: This figure depicts the values of the spectral index for all datapoints. They 
all lie on the same curve, with a dependence only on A, except for a small portion of the 
(g = 1) dataset. So most points really exist of 6 datapoints right on top of each other. In 
Ref. |20] there was a large dependence on the gauge parameter g, which now has totally 
disappeared. 



We assume that the cosmic strings are formed at the moment that the $ 
field leaves the inflationary valley. This happens if the process is relatively 
fast, in which case there will be a string production of approximately one 
length of string per Hubble volume. Figure (|19p shows at which angle the 
<3?-field will move from the S-axis. This angle always goes up very fast at 
the moment the fields leave the valley and then gradually goes down again. 
Since the leaving-angle is generally quite big, it is reasonable to assume that 
the strings are produced at the moment of going off-valley, but it is possible 
that the quantum fluctuations are still big enough to pull the field over 
the top of the Mexican hat for some time (9 . which would mean that the 
strings are effectively produced at some later instant, when the quantum- 
fluctuations become small as compared to the height of the Mexican hat. It 
would be interesting to do a precise analysis or simulation of this process. 

Normally one assumes that one length of cosmic string is produced per 
Hubble volume at S c . If there is inflation after the formation of strings the 
energy density stored in these cosmic strings will go down with the second 
power of the amount of inflatiorj^l. So, for example 30 off-valley efolds, 
would reduce the energy density of the strings with a factor e 60 ~ 10 26 . 

This certainly seems enough to get the string contribution below the 

14 Not the third power, because a cosmic string is a one dimensional object with a 
constant tension, which gets stretched by the expanding universe. 
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Figure 19: This figure denotes the maximum angle with respect to the S-axis of the steep- 
est descent path, for all datapoints. Generally the fields leave the S-axis quite abruptly 
and after that the angle gradually becomes smaller. 



10% bound, but this is not quite true. What matters is the string density 
at the moment of reionization. In between inflation and reionization the 
strings will interact and the string density will decrease. If the energy lost 
in this interaction goes into gravitational waves, then it does not contribute 
significantly to the WMAP anisotropy. As shown in Ref. [22] after a relevant 
cosmic string scale enters the Hubble radius again, it will go into a scaling 
regime, which makes the density of cosmic strings on a scale smaller than the 
Hubble radius approximately constant with respect to the Hubble radius. 
This means that the inter string distance grows with the Hubble radius. 
This is why the 10% string contribution bound only depends on £ [T7J: The 
density at reionization is always the same, but £ determines the tension 
of the string and therefore it determines the total energy-density stored in 
cosmic strings at reionization. The formula for the tension of a D-term string 
is very simple because it is a "critical coupling" or BPS string: 

H = 2vr£ . (39) 

Now if the number of off-valley efolds exceeds 60, there will be no string 
contribution for sure, because there simply are no strings. Up to about 5 
efolds less, the string contribution will still be very small, because the inter- 
string distance will not have entered the Hubble horizon at decoupling. For 
even smaller off-valley efolds the string contribution will go to the value as 
calculated before, due to the scaling. How fast the string contribution will 
go to this value depends on how fast the scaling will take place. 
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5.5 Pushing the top of the Mexican hat down 



As shown in section (13.41) . for a restricted set of parameters it is possible 
to lower the top of the Mexican hat of the fields in the true vacuum state. 
This could lead to strings with less tension and since the effect seems to be 
the largest for the parameters with a nice spectral index of 0.983, it is worth 
looking at. 

Since in the origin (S = 0, $ = 0) many of the particles are massless, 
we can actually do the analysis analytically. Evaluating all the masses, by 
looking carefully at tables ([3]), (jl]) and ([7]), we see that only four scalar 
masses contribute: M 2 (0 + ( 12 )) = g 2 £, and M 2 (0_( 12 )) = — <? 2 £. The rest 
of the fields are massless. Inserting this in the Coleman- Weinberg formula 
gives: 
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(40) 

We should compare this with the tree-level potential in the origin, see equa- 
tion Q. This gives a relative 1-loop correction in the origin of: 
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For our datapoints we used A m 
get: 
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S c . If we insert this into equation (I41h we 



This will give a negative 1-loop contribution for A < \^2Mpi, so all dat- 
apoints will have a 1-loop correction in the origin which counteracts the 
tree-level potential. Figure (f20|) shows the relative size of these 1-loop cor- 
rections. 

As you can see the 1-loop-contribution in the origin is very small indeed 
for almost all parameter combinations, but it can go up to about a maximum 
of 20% for values of the gauge parameter g very close to one. Of course we 
don't know what the two-loop contribution is going to do, so we cannot 
make a definite judgement to what extent this effect is important. The 
question remains open, whether the effect of this can be big enough to 
lower the string contribution below the 10% bound for the datapoints at 
large A. In any case, the best thing to do in order to get an expression 
for the tension of the strings, including this 1-loop correction, is to make a 
numerical model, find the string-solution and integrate the energy-density to 
get a corrected tension. Of course it will come out smaller than in equation 
(f39j) . but probably not small enough to get below the 10% string contribution 
bound. 
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Figure 20: Relative size of the negative 1-loop corrections in the origin (S = 0, <f>_ = 0) 
for all datapoints (a), zoomed in on the relevant region (b). For almost all choices of g we 
get a negligible 1-loop contribution. Only values of g close to one can give a significant 
contribution. 
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6 Conclusion 



In this paper we have looked in great depth into the D-term hybrid inflation 
model. The key ingredient we have added is a description of the fields 
after reaching the critical inflaton value and during spontaneous symmetry 
breaking, when the fields leave the inflationary valley. To analyse this we 
need a formula for the 1-loop quantum corrections to the potential energy 
density, outside the inflationary valley. In calculating this, problems arise 
related to the non-equilibrium Higgs mechanism. We presented a good way 
out of these problems, but there still is room for a thorough theoretical 
analysis of this phenomenon. 

Using our method, the parameter bounds of the D-term hybrid inflation 
model, compatible with WMAP measurements, improve a lot. A nice and 
simple relation between A and £ and independent of g is presented, that will 
give the correct CMB density perturbation. All parameter combinations 
with A < 4.5 • 10~ 5 give a cosmic string contribution less than 10% and 
also for values of A approximately smaller than 0.001 there is the possibility 
of diluting the strings after formation, which could also give a string con- 
tribution less than 10%. Moreover, where beforehand only g < 0.03 was 
compatible, now all values of g are possible. Results for the spectral index 
have no dependence on g. It lies between n = 0.983 and n = 1.0. The lower 
values are attained for approximately A bigger than 0.001. So the regions 
best compatible with the spectral index and with the cosmic string bound 
are exclusive. 

For D-term inflation to be compatible with WMAP, either n = 1.0 must 
prove to be the right value (it is currently at a 2a deviation), or neglected 
dynamical effects must give a contribution to the spectral index, making 
it smaller than one for A < 0.001, or 1-loop effects which lower the top of 
the Mexican hat potential in the true vacuum must lower the string tension 
below the 10% bound. 

In the future, work should be done related to the precise length scale 
of string formation in this model, including quantum fluctuations. The 
scaling behaviour can be examined, in cases where there are a lot of off- 
valley e- foldings. And also the influence of insertion of the true equations of 
motion, instead of the steepest descent path, on the spectral index might be 
interesting. Of course one could try to calculate the tension of the resulting 
strings, which have a lowered top of the Mexican hat potential. We saw that 
we could not neglect 1-loop corrections, so maybe even two- or higher-loop 
corrections may prove to be significant. But the first thing to do is to extend 
the used method to F-term hybrid inflation. 

So there are still a lot of improvements that can be made, but the model 
as a whole looks quite promising, which is something that could not be said 
before. 
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A Contribution from massless particles 



For some field-values, there will be massless particles contributing to the 
1-loop-correction in ([5]) . For example on the inflationary valley the fields 
and (S\S2 — Sis\)/\S\ will be massless. For these massless loops there is an 
indeterminate expression in ([5]), but taking a simple limit reveals that these 
massless particles really don't contribute to the loop-correction: 

hm (M 2 ln(M/A 2 m )) = 0, (43) 
where M is the mass-squared. 

B Contribution from imaginary masses 

We have identified all field-dependent masses of the physical particles in 
the region ($+ = 0, 3>_ and S arbitrary). In plugging them into equation 
([5]) there is still a problem. As can be seen in table j2j) the (mass) 2 -values 
are not always positive for all field values. Negative (moss) 2 -values plugged 
into the natural logarithm will give an imaginary contribution to the effective 
potential. This problem is resolved in Ref. |19j : 

Quantum-mechanically the effective potential can be interpreted as the 
energy-density of the quantum-state that minimises this energy-density, sub- 
ject to the condition that the field expectation values (in this case < <1>± > 
, < S >) are as given. So for every point in field-space you can define an ef- 
fective potential value. For classical potentials that are convex everywhere, 
these quantum-states will be homogeneous, meaning: concentrated around 
a point. But for partially concave classical potentials (which can give a 
spontaneous symmetry breaking) these quantum-states will be mixed states 
(because, if V" is negative, the mixed state will certainly have a lower energy 
than a homogeneous state). 

For our purposes we don't want to know the energy density of these 
mixed states, but rather we want to know the energy density of a homoge- 
neous state centred around the same concave point in field-space. Weinberg 
and Wu show in Ref. [19] that the encountered imaginary 1-loop-corrections 
to the effective potential correspond to these concave points (for mixed 
states). The nice thing is that the real part of the 1-loop correction can 
be interpreted as the 1-loop correction to the energy density of the homo- 
geneous stateEE For our calculations we should therefore only take the real 
part of the 1-loop effective potential. 

15 The imaginary part of the 1-loop corrections can be interpreted as a decay rate. 
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C Contribution from imaginary mass-squared val- 
ues 



As can be seen from tables ([2]) and ([!]) there is a possibility that some bosonic 
and fermionic (mass) 2 values are imaginary. This happens when the value 
under the square-root sign becomes negative. Therefore we'd better hope 
that this never happens. 

For the fermionic contribution it is easy to see that, indeed, E 2 -8\ 2 g 2 $ i 
is strictly positive for all field values, since it can be written as a sum of 
squares. 

E 2 - 8AV$ 4 = A 4 S 4 + 2A 2 (A 2 + 2g 2 ) S 2 $ 2 + (A 2 - 2g 2 f $ 4 > , (44) 

where S is shorthand notation for S 2 + S 2 . 

For the scalar contribution this is less trivial, because it is not clear how 
to write B 2 + AC as a sum of squares. 

B 2 + AC = {^-\ 2 S 2 - g 2 ^j + Q(A 2 - 3 5 2 )<D 2 + g 2 ^j 

+ ^ + \\ 2 9 2 )s 2 <$> 2 -g*e (45) 

However it is still possible that this function is strictly positive, because 
the three squares can not all be zero at the same time (fixing S and <3? can 
make two squares vanish, but then the third will not vanish). Calculating 
the minimum of this function in the (S, ^l-plane for a given A, g and £, this 
indeed turns out to be strictly positive 16 !. 

(^ 2 + ^) mm = ^ 2 (^)>0 (46) 

We conclude that there will never be imaginary mass-squared values in our 
computation, so we need not worry about this. 



D Higgs mechanism for polar and Cartesian fields 

Let us take a closer look at the definition of field-dependent particle masses 
in the case of a U(l) symmetry. Take the simplest possible global U(l) 
invariant Lagrangian density: 

JSf = (d^d^(j)-m 2 \(j)\ 2 . (47) 

16 The minimum is attained for: $ 2 = § jt^t and S 2 = - |] 
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We will look at it from two different coordinate-systems at exactly the same 
point in field space <fi = l/\/2$. The Cartesian coordinates 



4>= -J=($ + 0i + i0 2 ) (48) 
give the Lagrangian density 

2 = \%<fo&<h + \d^2^(t>2 - \m 2 ^\ + 4 + 2^i$ + $ 2 ) . (49) 

Both scalar fields have normal kinetic terms and both have a classical mass 
m, because the tree- level propagator (being the inverse of the quadratic part 
of the Lagrangian) in momentum space has a pole at k 2 = m 2 . If we look 
at the same Lagrangian in polar coordinates 

= -^(* + p)e* (50) 



we get 



1 — - 1 ,2/2 , o„^ , ^2n 



The />field enters in exactly the same way as the <^i-field before, but we can't 
read of the #-mass directly, because of the non-standard kinetic term, so 
this looks like a bad choice of coordinates. There is no term proportional to 
9 2 , because the second derivative of the potential along the circle trajectory 
corresponding to the #-field is zero. The real second derivative in the straight 
^-direction corresponds to the 02 _rriass an d is certainly not zero, see figure 
(|21|) . The trace of the (mass) 2 -matrix is just equal to the Laplaciarl^l. which 
is of course coordinate independent, as long as the coordinate trajectories 
are straight lines. This all seems obvious, and people will not interpret the 
#-field as a massless field in this case, because of the non-standard kinetic 
term. 

However if we make the global U(l) symmetry local, the gaugefield will 
enter. It's coupling with the #-field will cancel all ^-dependence and the 
masses can be easily evaluated again. Indeed using = + \d^6 the 
Lagrangian density becomes: 

2 = \dnpdPp + l -g 2 {<5> + pfB^ - \m 2 {<§> + pf - \f, u F^ . (52) 

The #-mass has disappeared. 

If instead we stay in the Cartesian basis, in which the masses were easy 
to read of before making the global U(l) invariance local, then the second 



17 The Laplacian is the sum of all second derivatives in perpendicular directions. 
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Figure 21: This picture shows a U(l) invariant potential. At the point A the second 
derivative following the #-circle-trajectory will be zero, exactly because of the U(l) invari- 
ance of the potential, but the genuine second derivative given by the intersection with a 
vertical plane will not be zero. Therefore the particle corresponding to the excitation of 
the field in this field-direction will not be massless. 



order term for the ' Goldstone f^l boson does not disappear. We can either 
take exactly the same definition of as before, then the Lagrangian density 
becomes* 

a> = l (^l(^ + ^l) + ^20l) 2 
2 (* + 0i) 2 + $ 

+ \g\{$> + <M 2 + 4)B^ - l -m 2 {^ + <M 2 + 4) - \F, U F^ . (53) 

Or we can use a locally equivalent definition of the Higgs-mixing = 
C » ~ |i 9 M^2 yielding: 

% = ^i^i+i 5 2 <I> 2 C^ , (54) 

where the dots stand for third- and higher order terms in the fields. In both 
cases the second order term for the Goldstone boson does not disappear in 
applying the Higgs mechanism. However we cannot associate a pole in the 
propagator with this term, because there is no kinetic term for the 02-field. 
Where in the global case one could either define the classical mass to be 
used in the Coleman Weinberg equation ([6]) by the pole in the propagator 

18 Maybe it is not right to call this a Goldstone boson, because this terminology is 
normally used in the global minimum (of, for example, the Mexican hat potential). In 
this minimum the Goldstone boson will always be massless, even for the global case. 
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or by the second derivative in the diagonal field direction. For the local 
case, including the Higgs mixing, both definitions no longer agree. In fact 
the second derivative in the 'Goldstone' field direction is the same as in the 
global case, but the 'Goldstone' mass coming from poles in the propagator 
arising from equation (|52H has vanished. 
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